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Abstract. The onset of surface fluidization of granular material in a vertically vibrated container, 
z = A cos (cut), is studied experimentally. Recently, for a column of spheres it has been theoreti- 
cally found Q that the particles lose contact if a certain condition for the acceleration amplitude 
z = Aco'^/g = f{Lo) holds. This result is in disagreement with other findings where the criterion 
z = Zcrit ~ const, was found to be the criterion of fluidization. We show that for a column of spheres a 
critical acceleration is not a proper criterion for fluidization and compare the results with theory. 
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1 introduction 

Granular material confined in a vertically vibrating con- 
tainer reveals cornplex effects such as surface structure 
formation (e.g. PP]), spontaneous heap formation (e.g. 
0), convection (e.g. |^) and others. A precondition com- 
mon to all these effects is that (at least) the particles at 
the free surface of the granular material lose contact to 
their neighbours for at least a small part of the oscillation 
period. If a granular material is agitated in a way so that 
the particles at the surface separate from their neighbours 
we will call the material fluidized. 

The conditions under which surface fluidization occurs 
are not clear yet. For sinusoidal vertical excitation of the 
container most of the literature reports that surface flu- 
idization starts as soon as the acceleration amplitude Auj'^ 
of the oscillation z = A cos tot exceeds gravity g or an- 
other critical constant. The Froude number is defined as 
r = AitP' I g and, hence, it has been reported in numer- 
ous publications that surface fluidization or, respectively, 
effects which require fluidization, occur for 7^ > 1. 

Few references report, however, that in numerical 
simulations surface fluidization was observed for F < 1. So 
far these results have not been confirmed experimentally. 

Whereas a rigid solid body, e.g. a single sphere on a 
oscillating surface, would certainly start to jump if the 
acceleration amplitude Auj^ exceeds gravity there are 
arguments which may lead to a different conclusion for 
a shaken amount of granular material, i.e., a many body 
system: 

— The deformation- force-law of contacting bodies may be 
nonlinear due to geometrical effects, even if the mate- 



rial deformation is small enough to assume linear ma- 
terial properties, i.e., Hooke's law. For contacting ide- 
ally elastic spheres it has been derived by Hertz W that 
the interaction force F scales with the deformation ^ 
as ^ ^1"^ . This law applies for all particle contacts 
(under mild conditions) provided the deformation ^ is 
small enough As a result, nonlinear phenomena like 
the propagation of solitary waves has been shown to 
appear in bead chains with Hertz contact law 
Under the influence of gravity the particles are de- 
formed differently even at rest, depending on their ver- 
tical position in the material. Therefore, the speed of 
sound in the material is not uniform but a function 
of the vertical coordinate. This property may lead to 
complicated pulse motion through the material and 
has been reported also for chains of contacting spheres 

0- . 

In a polydisperse granular material the geometrical 
properties, i.e., the contact network, varies with the 
applied force. If the material is loaded the particles are 
deformed and more and more contacts emerge. This ef- 
fect leads to a more complicated deformation - force 
law, F - with 1.5 < 7 < 4 [O. 



In a recent theoretical paper Q the granular material 
was modeled as a vertically shaken column of viscoelastic 
spheres. The main result of ||l| is that the column may be 
fluidized even if the condition F > \ does not hold. Instead 
a different condition for fluidization was derived: Assum- 
ing viscoelastic material properties the force between two 
contacting spheres of radius i? at vertical positions z-^ and 
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Zk+i reads 

Fk,k+i = -VR (MCfeffc+i + a^fc,fc+iA/Cfc,fc+i 



(1) 



where S,k,k+i = 2i?— \zk — Zk+i\ is the compression and ^ 
and a are elastic and dissipative material constants (de- 
tails see [|l2|). The height of the column of N spheres is 
L = 2NR. It was found that the top sphere separates from 
its neighbour if 
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where J-2/5 is the (complex) Bessel function and M ab- 
breviates 



M 
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(3) 



5 \ 9 J -\/3/i — 2iuja 
1 and g and p being gravity and material 



with i = 
density. 

For small frequency uj the Taylor expansion of the rhs 
of the condition (H) yields 



> 1 - B2UJ^ + Bauo^ 
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Both coefficients -82 and -84 are positive definite values, 
i.e., for all materials there is a range of frequency where 
the rhs. of the inequality (||) is smaller than one. This 
means that there exists always a frequency interval where 
the column fluidizes although the acceleration amplitude 
of shaking is smaller than gravity, Au}"^ / g < 1. The full 
derivation of the sketched theory and the discussion of 
the frequency range where the condition (^) holds can be 
found in Ref. For our experimental system consisting 
of steel spheres, of density p=7700 Kg/m'^, Young modulus 
Y—21Q GPa, Poisson ratio i/ ~ 0.29, in a column of L —0.6 
m, the Taylor expansion Eq. (^) holds for frequencies uj 
much smaller than 1000 s~^, which covers entirely the 
experimentally accessible range of frequencies. 

In the present paper we want to report experimental 
results on the onset of fluidization of a vertically shaken 
column of spheres, i.e., on the same system which was 
studied theoretically in Ref. . The dynamical behaviour 
of a vibrated one dimensional system of spheres in the well 
fluidized regime was studied e.g. injO,|l^. This regime is 
explicitely not considered here. In |15| | a column of parti- 
cles which interact as linear springs with linear damping, 
+ 2/i^ + i?^^ — 0, was studied with Molecular Dynam- 
ics. For their sets of parameters the authors found that 



fluidization occurs above F w 1.1, however, the authors 
claim that the onset of fluidization depends on the choice 
of the parameters f2 and p. Within their model it might 
be possible to get also fluidization for _r < 1, depend- 
ing on the material parameters. If one considers spheres, 
however, in the limit of vanishing damping the Hertz law 
implies that the duration of contact behaves as tc ^ 9~^^^, 
whereas for the model considered in [p^ one finds a con- 
stant tc ~ const. = tt/I?. Hence, the linear spring model 
certainly does not describe the contact of spheres, not even 
in the case of pure elastic contact. For the same linear 
spring system in [l^ the detachment effect was reported, 
the same effect for viscoelastic spheres was discussed in 



2 Experimental setup 

We expect that the effect to be measured is quite small, 
i.e., for the critical parameters of shaking when the mate- 
rial starts to fluidize we expect 1 — AuP' / g to be a small 
value. Therefore, we have to adjust the amplitude of shak- 
ing A and in particular the frequency w with good accu- 
racy. Moreover, we have to assure that no energy originat- 
ing from different frequencies than the driving frequency 
enters the column of spheres, e.g., due to high frequency 
noise. 

The column of spheres was confined in a framework of 
three aluminum bars, which assures pure one dimensional 
motion of the spheres at low friction, see Fig. |l|. 



view from above 



front view 




Fig. 1. The column of spheres moves in a framework of 
aluminum bars 



The framework with the spheres was mounted on a pre- 
cisely vertically balanced linear bearing which was driven 
by a stepping motor via a crankshaft with adjustable ec- 
centricity. Fig. ^. 

The motor was computer controlled, i.e., in precisely 
20,000 impulses the motor axis revolves once. This high 
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angular resolution provides quasi-steady motion. The mo- 
tion was smoothcned additionally using a flywheel (with 
diameter 160 mm, thickness 32 mm, mass 4.5 kg and mo- 
ment of inertia J — 0.0144 kgm^) which was fixed on the 
axis and a hard rubber clutch connecting the motor axis 
and the axis of the shaking device (see Fig. ||). Using an 
acceleration sensor we checked that the finite step size per 
computer signal does not influence the sinusoidal motion 
of the column of spheres. The entire mechanical device was 
mounted on a massive oscillation damping table. Figure ^ 
sketches the device and Fig. || shows a photograph. 

In the experiments the column consists of up to 20 ball 
bearing balls of radius R = 12.5 mm. There are several 
possibilities to determine the onset of fluidization, i.e., for 
fixed amplitude A the critical frequency to when the top 
sphere separates from its neighbour: 

— Direct observation: a comoving camera was mounted 
to the probe carrier and a LASER pointer was fixed at 
the opposite side of the column. As soon as fiuidization 
starts at the critical frequency the emerging submil- 
limeter gap between the top particle and its neighbour 
can be recognized as a flash on the screen. 

— Acoustic method: When fluidization occurs there emer- 
ges a periodic acoustic signal originating from the col- 
lision of the top particles. This signal can be detected 
either using a microphone or the acceleration sensor 
which was mounted on the top of the device, see Fig. ||. 

— Electric circuit: We attached a very thin wire (diam- 
eter < 0.1mm) to the top particle and its neighbour 
and observed the contact of the spheres using an oscil- 
loscope. 

We checked all three methods and found that the re- 
sults are very similar. We have chosen the electric method 




Fig. 4. Experimental device with column of steel spheres 
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Fig. 5. The onset of fluidization was determined by observ- 
ing interruptions of tlie electrical circuit. 



since it turned out to be the most simple and most reliable 
one. The direct observation requires complicated mechan- 
ical tuning and for the acoustic method one has to filter 
the periodic signal from other sources of sound as, e.g., 
the sound of the motor. 



3 Raw Data 

For fixed amplitude A we want to determine the accord- 
ing smallest frequency uj for which the top sphere sepa- 
rates from its neighbour in each period. Figure || shows a 
sketch of the oscilloscope signal for subcritical (top), crit- 
ical (middle) and over-critical frequency, all for the same 
amplitude of shaking. For subcritical frequency the top 
particle follows the motion of the vibrating table resulting 
in a sinusoidal signal of the acceleration sensor and a hor- 
izontal line for the conductivity (limited by a resistor). At 
over-critical frequency the conductivity drops to zero for 
a certain time interval in each period. Close to the critical 
frequency the electrical contact between the spheres is not 
in each period interrupted, instead we get a noisy signal. 

To collect data points we tuned the adjustable crank- 
shaft to a desired amplitude A. Then the motor was ac- 
celerated up to a rotation speed which is slightly less than 
the speed where we expect fluidization (initial ramp). At 
this value all spheres are permanently in contact (top in 
Fig. H). From this point on we increased the motor speed 
in very small steps until we observe separation of the top 
sphere in each period (bottom in Fig. 6|). This procedure 
was repeated 10 times for all amplitudes independently for 
2 and 20 spheres. The critical frequency for each ampli- 
tude was determined as the average over the 10 indepen- 
dent measurements. Figure |7| shows the raw data points 
with error bars. 

The somewhat unusual representation of the data as 
acceleration vs. amplitude is due to the fact that our ex- 
perimental input is the amplitude, the output, as stated 
above, is the according critical frequency and hence the 
acceleration. 



1 1 ; 


1^ 

1 1 1 1 1 1 


T-l 1- 

1 1 1 
( 1 1 

1 1 1 


1 r-i r-i r-i H- 

1 1 1 1 1 1 
1 1 1 I 1 1 1 
1 1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 

1 ' ', ' 1 ' 

1 1 1 


I 1 i 1 i 1 1 1 

II 1 1 1 1 1 1 
11 II 1 1 1 1 
II 1 1 1 1 1 1 
11 II 1 1 1 1 
II 1 I 1 1 1 1 
II 1 1 1 1 1 1 
11 1 1 1 1 1 1 
II 1 1 1 1 1 1 

V V V V \l 



time 



/\ l\ /^ /I /^ 

M I 1 M I 1 I 1 



I I 
1 I 
1 I 



I I 
1 I 
I I 



I I 1 
1 I I 



I 1 I 

-I 1 r 



1 I I I 1 

I I 1 I 

I I 1 I 

1 I I I 

I I 1 I 

1 I 1 I 

1 I 1 I 

I I 1 I 

time 



I I 
1 I 
I I 



I I 
I I 
w 

time 



' I 
' I 
' I 
' I 
' I 
' I 
' I 

1 I I 
1; 



Fig. 6. For subcritical sliaking the particles are always in con- 
tact, i.e., the conductivity (solid line) is constant (top). For 
over-critical shaking the conductivity drops to zero in each 
period (bottom). For critical values the conductivity shows ir- 
regular peaks (middle) 



In Figure || the critical acceleration of the container 
{Auj'^) vs. frequency uj for 2 and 20 beads is shown. One 
can see that even for two beads the critical acceleration, 
which is supposed to be a constant varies significantly with 
the frequency, i.e., it decreases with increasing frequency. 
The curve for 20 beads shows a similar behaviour, how- 
ever, the critical acceleration for 20 beads is smaller than 
the critical acceleration for 2 beads. This is a first indi- 
cation to an amplification effect similar as theoretically 
predicted. The experimental device reveals its own res- 
onances, namely 7.6 Hz for two beads and 11.8 Hz for 
twenty beads. These resonances are caused by the me- 
chanical properties of our experimental setup but not by 
the properties of the columns of beads themselves. Data 
points which belong to these frequencies have been elimi- 
nated for further analysis (see Fig. |^). 

The data points are affected by large statistical errors, 
i.e., the data points scatter around a smooth curve which 
one expects to find if only systematic effects were relevant. 
Therefore, we have to apply a sophisticated data analy- 
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Fig. 7. All measured data points. In the upper figure er- 
ror bars have been attached to the measurement for a single 
free moving bead (and another one serving as a well defined 
lower boundary). In the lower part the same data is shown 
with error bars for the measurement using 20 spheres. Data 
points marked by dashed lines have been disregarded for fur- 
ther analysis since they are adulterated by resonances of the 
equipment (see text). 
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Fig. 8. The critical acceleration versus frequency. Circles 
represent data of two spheres (reference), diamonds represent 
data for 20 spheres. 



sis to extract the amplification effect in order to compare 
it with the theoretical prediction. This procedure is de- 
scribed in the following sections. 



4 Discussion of the Measuring Procedure 

The measurement is affected by different systematic and 
random errors: 

— Due to the generation of the shaking motion via a con- 
necting rod of finite length {I = 18 cm) there is a 
systematic deviation of the motion from the ideal si- 
nusoidal shaking. Instead of a sinusoidal oscillation the 
generated oscillation is 



z = A cos + — sin^ ( 
For small ratios A/ 1 we can write 

A' . . 



I. 



z K, A cos ( 



2L 



sm 



(5) 



(6) 



This leads to a constant shift in the extremal acceler- 
ation of the container 



^Extr 



1j 



AJ' ±1 



We see that the absolute value of the acceleration at 
the upper turning point is too large hy A/l which can 
be as much as 3%. 

Uncertainty in determining the driving frequency and 
amplitude are very small as described above. 
Duo to the noise caused by the motor and the fric- 
tion of the linear bearings the systematic (ideally si- 
nusoidal) vertical motion of the column is superposed 
by a spectrum of high frequencies. These waves are 
source of an undesircd extra energy feed which is not 
controllable and causes errors. Measurements of this 
noise using the acceleration sensor have shown that it 
does not depend much on the height of the column, it 
is determined mainly by the motor and the bearings. 
This property suggests to measure our effect not di- 
rectly but to compare the results for a column of 20 
steel spheres with the results for a "column" of two 
spheres (see below). 



5 Data Analysis 

To eliminate systematic errors we compare the result for a 
column of 20 spheres with the equivalent measurement for 

a "column" of only two spheres. The column of two spheres 
serves as a reference. Actually, the second bead is fixed to 
the bottom and oscillates rigidly while the top bead is free. 
Effectively there is only one free moving sphere, for one 
single sphere taking off from a flat wall could not serve as 
a reference for the case when the top bead of a column of 
N spheres separates from the rest. This approach is based 
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on the assumption that the systematic errors originating 
from the motor noise and from the bearings affect both 
systems almost in the same way. With the above men- 
tioned method (see section ||) we determined the critical 
frequency of driving at which the topmost bead starts to 
jump for a given shaking amplitude. From this we cal- 
culated the critical acceleration. Due to the influence of 
systematic errors also the curve acrit('^) for two beads, 
which is supposed to be constant, shows a significant de- 
pendence on the driving frequency, resp. amplitude. As 
already mentioned before, in Fig. H one can see that the 
critical acceleration decreases with increasing frequency. 
The curve for 20 spheres shows a similar behaviour but 
the critical amplitude is always smaller than the critical 
amplitude for 2 spheres at the same driving frequency. To 
identify the amplification effect, i.e., to eliminate the error 
of the apparatus we divide the critical acceleration for 20 
spheres by the critical amplitude for 2 spheres at the same 
frequency which yields the absolute value for the critical 
Froude number for 20 spheres. This method is applicable 
if all significant systematic errors affect the measurements 
by a factor which is independent on the number of beads. 
This precondition holds certainly for the influence of the 
limited length of the connecting rod and holds to a rea- 
sonable degree for the influence of the higher harmonics 
in the driving motion due to the discrete steps in the mo- 
tor motion. Since we believe these errors to be the most 
significant ones our method should give reliable results. 

Following this idea we fit both sets of data, for 2 and 
20 spheres, to the function 



/(cj) = (020;^ -I- aiUJ + ao) 



(8) 



and take the ratio of these curves. There are other possi- 
bilities to construct fit-functions which reflect the desired 
properties of monotonous decay and significant curvature, 
it turns out that the final result does not depend sensi- 
tively on the choice of this function. 
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Fig. 9. All relevant data points together with a least squares 
fit curve according to template Eq. Full circles: 2 beads, 
diamonds: 20 beads. 



The ratio curve, which represents the predicted am- 
plification effect, can be described in good approximation 
as a decreasing parabola with maximum at w = 0. Its 
curvature is small, in agreement with the expansion given 
by Eq. (^. Dissipative effects, represented by the mate- 
rial constant a, are not observable at these frequencies, 
as the restitution is high for all velocities (0.90 or more): 
we estimated the second term of -B4 being of the order 
of 10^^ or less. This curve, multiplied by the reference 
acceleration 9.81 m/s^, is shown in Fig. ^ The ratio 
between the fitting curves in Fig. | is not exactly 1 at 
a; = 0, but by about 3% smaller. This discrepancy is due 
to the above mentioned systematic errors, namely due to 
the higher harmonics in the motion of the stepping mo- 
tor (see Fig. |ll|) which can have different magnitude at 
different load. Therefore the division of the curve for 20 
beads by the reference curve alone can not completely re- 
move this systematic error. Other systematic errors such 
as deviations from an ideal sinusoidal motion due to the 
finite length of the connecting rod have been removed by 
the data analysis procedure. 
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Fig. 10. Dashed line: the ratio between the fitted curve for 
20 and 2 beads (see Fig. ^), multiplied by the theoretical 
value 9.81 m/s^. Solid line: theoretical values of the ampli- 
fication factor given by Eq. (^) (upper curve) and the same 
result translated vertically for comparison with the experi- 
mental one (bottom solid curve) . The constant shift between 
both solid lines corresponds to about 3% of the vertical axis 
due to the systematic errors discussed in the text. The dot- 
dashed line is the acceleration of gravity, plotted for refer- 
ence. 



The influence of the higher harmonics can be observed 
in Fig. |l^ showing the output of the acceleration sensor on 
top of the column. The sinusoidal motion of the container 
with driving frequency 10 Hz is superimposed by an os- 
cillation of frequency of about 150 Hz. When calculating 
the maximum acceleration of this motion one gets a value 
which is up to 20% higher than AuP . Thus, by taking Alu"^ 
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as the maximum acceleration of the container we under- 
estimate the actual acceleration, which, since the magni- 
tude of the high frequency oscillation varies with different 
load, contributes to the difference between the maximum 
of the ratio curve in Fig. |l^ and the theoretical value of 
9.81m/s2. 

6 I 1 



4 - 




0.0 0.1 0.2 0.3 0.4 

time in s 



Fig. 11. The output of the acceleration sensor (lower curve). 
One can see that the sinusoidal curve of the driving frequency 
is disturbed by higher frequency oscillations. The higher fre- 
quencies account to up to 20% of the maximum acceleration. 
The upper curve shows the conductivity measurement of the 
topmost beads. The peaks correspond to jumping of the top- 
most bead. 



6 Discussion 

We investigated experimentally the onset of fluidization of 
a vertically vibrated column of spheres. It was found that 
one observes fluidization even if the amplitude of acceler- 
ation of the vibrating motion is smaller than the accelera- 
tion due to gravity, AuP' /g < 1. This result is in agreement 
with the theoretical prediction 0. The quantitative the- 
oretical result for the critical parameters for the onset of 
fluidization given by Eqs. (||) and (^) does not contradict 
the experiment but could also not be conclusively con- 
firmed by the experiment due to insufficiencies of the ex- 
perimental setup. The most important shortcoming is the 
existence of high frequency vibrations mainly due to the 
motor noise and the bearings. While qualitatively agree- 
ment is shown, due to these limitations the experimental 
measurements are not completely conclusive to verify the 
theoretical results: a quantitative check of the theoreti- 
cal predictions requires a more sophisticated experimen- 
tal setup, directed to reduce noise and enhance the effect. 
The noise can be reduced by improving mechanical isola- 
tion: enclosing the column of spheres into a massive block 
or wall, using very rigid bearings and separating the mo- 
tor at a good distance from the column. Higher columns 



or more balls will help to reduce the noise, but one has 
to work harder on a correct vertical alignment; for this 
reason the beads should not be very small. The use of a 
softer material alone would enhance the effect, but then 
one would like still good restitution, smooth surface and 
good conductivity (for the electrical method of determina- 
tion) , so the choice of material is not straightforward. Also 
higher frequency will enhance fluidization at lower Froude 
numbers, but then controlling amplitude (and noise) will 
become much more critical. 



We thank Hans Scholz for discussion and Christine Rosinska 
for technical aid. 



References 

1. T. Poschel, T. Schwager, and C. Saluefia, Phys. Rev. E 62, 
1361 (2000). 

2. C. A. de Coulomb, Memoires de Mathematiques et de 
Physique Presentes a L'Academie Royale des Sciences par 
Divers Savans et Lus dans les Assemblees 7, 343 (1773). 

3. F. Melo, P. Umbanhowar, and H. L. Swinney, Phys. Rev. 
Lett, bf 72, 172 (1994); T. Metcalfe, J. B. Knight, and H. 
M. Jaeger, Physica A 236, 202 (1997). 

4. F. Dinkelacker, A. Hiibler, and E. Liischer, Biol. Cybern. 
56, 51 (1987). 

5. J. B. Knight, E. E. Ehrichs, V. Yu. Kuperman, J. K. Flint, 
H. M. Jaeger, and S. R. Nagel, Science 267 1632 (1995). 

6. G. C. Barker and A. Mehta, Nature 364 486 (1993); J. A. 
C. Gallas, H.J. Herrmann, and S. Sokolowski, Physica A 
189, 437 (1992); Y-h Taguchi, Phys. Rev. Lett. 69, 1367 
(1992). 

7. H. Hertz, J. f. reine u. angewandte Math. 92, 156 (1882). 

8. N. V. Brilliantov, in preparation. 

9. C. Coste, E. Falcon, and S. Fauve, Phys. Rev. E 56, 6104 
(1997). 

10. R. S. Sinkovits, and S. Sen, Phys. Rev. Lett, bf 74, 2686 
(1995); S. Sen, M. Manciu, and J. D. Wright, Phys. Rev. 
E 57, 2386 (1998). 

11. D. Stauffer, H. J. Herrmann, and S. Roux, Europ. Phys. 
J. 3, 265 (1987). 

12. N. V. Brilliantov, F. Spahn, J.-M. Hertzsch, and T. 
Poschel, Phys. Rev. E 53, 5382 (1996). 

13. B. Bernu, F. Delyon, and R. Mazighi, Phys. Rev. E 50, 
4551 (1994). 

14. S. Luding, E. Clement, A. Blumen, J. Rajchenbach, and 
J. Duran, Phys. Rev. E 50, 1634. 

15. ibid., p. R1762. 

16. ibid., p. 4113. 

17. S. Luding, PhD-Thesis (1994). 



